Intensity fluctuations in bimodal micropillar lasers enhanced by quantum-dot gain 
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We investigate correlations between orthogonally polarized cavity modes of a bimodal micropillar 
laser with a single layer of self-assembled quantum dots in the active region. While one emission mode 
of the microlaser demonstrates a characteristic s-shaped input-output curve, the output intensity 
of the second mode saturates and even decreases with increasing injection current above threshold. 
Measuring the photon auto-correlation function g^-^^r) of the light emission confirms the onset of 
lasing in the first mode with g^^\0) approaching unity above threshold. In contrast, strong photon 
bunching associated with super-thermal values of g^'^\0) is detected for the other mode for currents 
above threshold. This behavior is attributed to gain competition of the two modes induced by the 
common gain material, which is confirmed by photon crosscorrelation measurements revealing a clear 
anti-correlation between emission events of the two modes. The experimental studies are in excellent 
qualitative agreement with theoretical studies based on a microscopic semiconductor theory, which 
we extend to the case of two modes interacting with the common gain medium. Moreover, we 
treat the problem by an extended birth-death model for two interacting modes, which reveals, that 
the photon probability distribution of each mode has a double peak structure, indicating switching 
behavior of the modes for the pump rates around threshold. 



I. INTRODUCTION 

Quantum dot - microcavities are a very attractive sys- 
tem to study quantum optical effects in the solid state^. 
Apart from research on fundamental light matter inter- 
action in the weak and strong coupling regime of cavity 
quantum electrodynamics^^^, they offer the possibility 
to investigate stimulated emission in a regime approach- 
ing the ultimate limit of a thresholdless laser based on 
a single zero-dimensional gain center—. Studies in this 
field include, e.g., technological works on optically and 
electrically pumped microlasers aiming at an increase of 
the /3-factor which expresses the fraction of spontaneous 
emission coupled into the lasing mode^"— . In high /3- 
microlasers it becomes increasingly difficult to identify 
the transition from spontaneous emission to stimulated 
emission at threshold via their input-output character- 
istics". This issue has triggered comprehensive experi- 
mental and theoretical research activities on the photon 
statistics of emission in terms of intensity autocorrela- 
tion function in order to unambiguously identify the on- 
set of stimulated emission at thresholdSiH^— . Moreover, 
the autocorrelation function is very beneficial to iden- 
tify single quantum dot controlled lasing eSecta^'^^'^ and 
to reveal other intriguing effects such as correlations be- 
tween individual photon emission events^ and chaotic 
behaviour of feedback coupled microlasersiS. 

The research efforts on microcavity lasers so far have 
focused mostly on emission features based on the inter- 
action between a single laser mode and the quantum dot 
gain medium. Going beyond this standard investigations, 
micropillar lasers with a bimodal emission spectrum al- 



low one to address the coupling of two orthogonal optical 
modes via the common gain medium which can lead to 
characteristic oscillations in the coherence properties'^ 
and an enhanced sensitivity on external perturbations in 
the presence of optical self- feedback-^ . 

In the present work, we perform a detailed experimen- 
tal and theoretical analysis of the mode coupling and gain 
competition of bimodal, electrically pumped micropillar 
lasers. A convenient measure for the study of the sta- 
tistical properties of the electromagnetic field emission is 
the set of intensity correlation functions: 
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(1) 



where — 1j2, with delay time r and photon annihi- 
lation operators bi and 62 of the mode 1 and the mode 
2, correspondingly. The gain competition is reflected in 

distinct differences in the input-output characteristic and 

(2) 

the autocorrelation function g^^ (r) of the two optical 

(2) 

modes. Moreover, the crosscorrelation function §{2 (t) 
can illustrate correlations between emission events from 
the two modes. In order to describe and analysis these 
specific features of bimodal microlasers we extend the 
microscopic semiconductor modeli^ accordingly by tak- 
ing mode interactions into account. Similarly, we extend 
a standard birth-death approach^^ for the description of 
bimodal lasers. While the microscopic semiconductor 
theory is applied to model the input-output character- 
istics, the intensity correlation functions of the laser and 
the gain competition between the two emission modes 
within a strict mathematical framework, the extended 
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birth-death approach aUows for a more intuitive under- 
standing of the underlying photon statistics. 

The paper is organized as fohows. In section HIl the ex- 
perimental results obtained from an electrically pumped, 
bimodal micropillar laser will be presented. Section IIIII 
deals with the theoretical description of the experimen- 
tal data and is divided into two subsections addressing 
a microscopic semiconductor theory, and an extended 
birth-death approach, respectively. The paper closes in 
section IIVI with a comparison of the experimental and 
theoretical results and a conclusion. 



II. EXPERIMENT 

The electrically pumped micropillar lasers are based 
on planar AlAs/GaAs microcavity structures which in- 
cludes an active layer consisting of a single layer of 
Ino.3Gao.3As. High resolution electron beam lithogra- 
phy, plasma enhanced etching and metal deposition have 
been applied to fabricate high quality electrically pumped 
microlasers. For more details on the sample processing 
we refer to Ref.^^. The microlasers have been investi- 
gated at low temperature (20 K) using a high resolution 
micro-electroluminescence (/iEL) setup. A linear polar- 
izer in combination with a A/4- wave-plate is installed 
in front of the entrance slit of the monochromator in 
order to perform polarization resolved measurements of 
the laser signal. The photon statistics of the emitted 
light has been studied by means of the measurement 
of the photon autocorrelation function g^^^'{T), that has 
been carried out using a fiber coupled Hanbury-Brown 
and Twiss (HBT) configuration with a temporal resolu- 
tion Tirt — 40 ps. The HBT configuration is coupled to 
the output slit of the monochromator which has a focal 
length of / = 0.75 m. The interaction of the orthogo- 
nally polarized modes of the microlaser has been investi- 
gated by means of photon crosscorrelation measurements. 
For this purpose, the light emitted by the microlaser is 
split by a polarization-maintaining 50/50 beamsplitter 
and coupled into two monochromators (/ — 0.75 m), each 
of which is equipped with a linear polarizer at the input 
slit and a fiber coupled single photon counting module at 
the output slit. This configuration allows us to perform 
polarization resolved crosscorrelation measurements with 
a spectral resolution of 25 ^eV. 

First, let us focus on the input-output characteristics 
of the microlaser. Due to slight asymmetry of the cross- 
section of the pillar and the ring-shaped contact the de- 
generacy of the fundamental mode in the pillar micro- 
cavity is lifted and two distinct linearly polarized modes 
are supported--. In this context, the spectral splitting 
A12 and accordingly the overlap between the two modes 
plays an important role for the studies of emission of bi- 
modal cavities. Figure [T] shows representative polariza- 
tion resolved spectra of an electrically pumped bimodal 
microlaser at threshold (injection current, linj =5.1 fJ-A). 
The two linearly polarized modes are split in energy by 
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FIG. 1. (color online) Polarization resolved fiEL emission 
spectra of a microlaser with a diameter of 3 fim. The elec- 
tromagnetic field emission features two orthogonally polar- 
ized cavity modes, the mode 1 (Q = 13900) and the mode 2 
(Q — 13100) with a spectral separation of 103 fieV (Injection 
current: Ii„j = 5.1 /^A). 

103 //eV and have absorption limited Q-factors oi Q — 
13900 (mode 1) and Q = 13100 (mode 2) at the thresh- 
old. The input-output characteristic of the bimodal mi- 
crolaser is presented in Fig.[2ja). We observe pronounced 
differences between the two modes: while mode 1 shows 
a standard "s" -shaped input-output characteristic with 
a threshold current of about Ith = 5.1 /lA, the intensity 
of mode 2 saturates at /Ith = 2 and even drops down 
for injection currents exceeding I^j / Ith = '^■5- This be- 
havior indicates a pronounced competition between the 
modes 1 and 2 which is mediated by the common QD gain 
material as it will be further elaborated in the following. 

Further, to study the lasing features we extract the 
emission linewidths of the two modes and plot them as 
a function of the injection current in Fig. ^h). The 
linewidths of the modes 1 and 2 have similar magni- 
tude and decrease strongly at threshold which reflects 
enhanced temporal coherence in the lasing regime. In- 
terestingly, while the linewidth of the mode 1 stays at 
a resolution limited value of 25 fJ,eV, a slight increase 
of the linewidth can be observed for the mode 2 above 
linj /Ith = 3. This is in agreement with the decreasing 
emission intensity seen in Fig. [2{a) , which indicates an 
increasing contribution of spontaneous emission in the 
mode 2 at high injection currents. 

In order to verify the interpretation of mode coupling 
in terms of gain competition, we have performed cross- 
correlation measurements between the modes 1 and 2 at 
different injection currents. An illustration of such a mea- 
surement is presented in Fig. [2l^e) for Imj/hh = 3. The 

(2) 

cross-correlation function (t) shows a pronounced dip 

(2) 

912 min — 0-62 at T = which indicates an anti-correlation 
between emission events from the two laser modes. 
The anti-correlated emission occurs at a characteristic 
timescale of ti2 = 3.8 ns. Figure [2l^d) reveals that the 
crosscorrelation function 5^2 (0) strongly depends on the 
injection current. It is useful to note, that, as it is seen 
from Fig. [2] in the regime of certain injection currents 
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FIG. 2. (color online) Experimental characteristics of a bi- 
modal micropillar laser with a diameter of 3 fim. (a) Input- 
output characteristic, (b) emission mode linewidth and the 
photon (c) auto- and (d) crosscorrelation functions g\^{0), 
g^{0) and (0) of emission from modes 1 and 2, respec- 
tively. Panels (e) and (f) show exemplary crosscorrelation 
g['^{T) and autocorrelation g^ (t) measurements at an injec- 
tion current of 



above the threshold {2.1 < Imj / Ith < 3.3), the intensity 
of mode 2 decreases and the statistics of mode 2 demon- 
strates strongly super-Poissonian behavior, whereas the 
anti-correlation between the modes is the strongest. 

The interplay between the two emission modes is also 
accompanied by strong temporal intensity fluctuations 
which are identified by measuring the photon autocor- 
relation function of the two competing modes for differ- 
ent injection currents. The respective dependencies, i.e. 
Qiii^) and 522 (0) versus injection current, are plotted 



in Fig. [U^c), while Fig. [2jf) shows the autocorrelation 

(2) 

function 522 (''') ^^r Imj/Ith = 3. The mode 1 shows the 

(2) 

typical maximum of g-n (0) around threshold, which in- 
dicates the transition from spontaneous emission to stim- 

(2) 

ulated emission, where 5ii (0) is lower than expected 
from theory due to the limited temporal resolution of the 
HBT^**. In contrast, the autocorrelation function of the 
mode 2, 522 (0) increases strongly at the pump rates well 
above threshold and reaches a maximum value of 3.08 

at hnj/Ith — 3. This value is significantly higher than 

(2) 

522 (0) — 2, expected for thermal light and, therefore, 
can not be explained by standard photon statistics. 

It is important to note, that similar statistical proper- 
ties of the emission, i.e. strong super-Poissonian behavior 
for the weak mode, has been also observed for microlasers 
in the presence of an external mirror, where a delayed 
feedback of the emitted signal disturbs laser operation 
and leads to strong bunching for the weak modei^. 



III. THEORY 

To develop a theoretical framework for the study of 
the coupled carrier-photon system in the bimodal cavity 
we consider two different theoretical approaches. In the 
first, a microscopic theory of light-matter interaction of 
semiconductor QDs with the cavity field is given, which 
allows the derivation of the equations of motion for quan- 
tities of interest. In the second approach, starting with 
the master equation, statistics of the photon distribution 
can be derived for the case of two-level carriers. 



A. Microscopic Semiconductor Theory 

To study the interaction of QDs with the electromag- 
netic field inside an optical bimodal microcavities we have 
extended the Microscopic Semiconductor Theorj*i^ to the 
case of two modes and photon crosscorrelation functions. 



1. Physical Model 

The Microscopic Semiconductor Theory allows for in- 
clusion of many-body effects of the carriers and can be 
used to calculate correlations required to determine the 
statistics of the emission of microcavities with active QDs 
(for a review see, e. g., Ref J^). The calculations are based 
on the cluster expansion truncation scheme of the equa- 
tions of motion for operator expectation values^^. 

In what follows we assume that only two confined QD 
shells for both electrons and holes are relevant: whereas 
the resonant interaction with the electromagnetic field of 
the bimodal cavity is due to the coupling with the s-shell 
transition, the carrier generation due to electrical pump- 
ing is to take place in the p-shell. The assumption suits 
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well also for an experimental situation, where the elec- 
trical pumping is to take place via injection of electrons 
and holes into the wetting layer and subsequent fast re- 
laxation to the discrete electronic states of the QDs. Fur- 
ther, carrier-carrier and carrier-phonon scattering contri- 
butions to the dynamics are evaluated using a relaxation 
time approximation, where the relaxation towards quasi- 
equilibrium is given in terms of a relaxation rate^^. 

To be more specific, let us consider a bimodal micro- 
cavity with the QDs as gain medium with the driving 
performed by the recombination of carriers in the valence 
and conduction bands. The Hamiltonian that governs the 
temporal evolution of the overall system can be given in 
the form 



(2) 



where -ffcaiT the single-particle contributions for con- 
duction and valence band carriers with the energies e^'", 

3 J 

and the two-particle Coulomb interaction is given by2i 



k'jj'k 



E 

k'jj'k 



yk^jfkcWjVj'Ck- 



(4) 



In the above, Cj (cj) and Vj (wj) are fermionic opera- 
tors that annihilate (create) a conduction-band carrier 
in the state \j)c and a valence-band carrier in the state 
\i)v, respectively. Further, the Hamiltonian of the elec- 
tromagnetic field modes inside the cavity reads 



(5) 



where is the bosonic annihilation (creation) oper- 

ator of the ^th mode of the cavity. 

The energy of interaction of the QDs with the electro- 
magnetic field inside the cavity in dipole approximation 
can be given by: 



jb^) + }l.c. 



(6) 



where the approximation of equal wave-function en- 
velopes for conduction- and valence-band states is used. 
Moreover, for simplicity the coupling strength g^j is as- 
sumed to be real. 

The Hamiltonian given by Eq. ^ together with 
Eqs. ([3]-[6]) determines the dynamical evolution of the 
carrier and field operators and, in particular, the time 
evolution for operator expectation values. The equations 
of motion for quantities of interest, as for example the 



average photon number in the cavity modes and the av- 
erage electron population in the conduction and valence 
bands, have source terms that contain operator expec- 
tation values of higher order. In this way, the approach 
bears an infinite hierarchy of equations of motion for var- 
ious expectation values for photon and carrier operators. 
To perform a consistent truncation of the equations the 
cluster expansion scheme is applied (for details, see RefJ^ 
and references therein). Namely, starting from the ex- 
pectation values of the first order of photon operators, 
the equations of motion for operator expectation values 
are replaced by equations of motion for correlation func- 
tions. For example, instead of the equations of motion for 
expectation values of amplitudes of the cavity mode op- 
erators the equations of motion for corresponding 

amplitude correlation functions — (^^^c) ~(^|)(^c) 

are used. Then, to achieve a consistent classification and 
inclusion of correlations up to a certain order the trun- 
cation of the equations for correlation functions rather 
than for expectation values is performed. 

In particular, in the case of a system without coherent 
external excitation = = and {c^Vj') = hold. 
Therefore, applying rotating- wave approximation here 
and thereafter, Heisenberg equations of motion for am- 
plitude correlation functions of the mode operators can 
be given by 



dt 



Siblbc) 



T.{9iA'^hh)+9eAv]c4)), (7) 



where is the loss rate of the ^th cavity mode and 
q= 1 . . . N, with N being the total number of QDs. 
Note, that both cavity-mode amplitude correlation func- 
tions 5{b^^b^) and the coupled photon-assisted polarization 

amplitude correlations d{vjcjb^^) and 6{CjVjbQ) are classi- 
fied as doublet terms in the cluster expansion scheme, 
i.e., they correspond to an excitation of two electrons 
(four carrier operators). The equation of motion for the 
photon-assisted polarization amplitude correlation read 
[see also Eq. (|Al|) in Appendix lA] : 

^Siv^ajbl) = -iiA^j - in^ ~ ir)S{v]cjbl) 



+9i'AcyA'k^ - 9i'A4^A,b^)] ' (8) 



where A^j = i{ej — Ej) — hjj(_ is the detuning of the 
^th cavity-mode from the QD transition and F is a phe- 
nomenological dephasing parameter describing spectral 
line broadening. In the case of a bimodal cavity only 
the cavity modes with indices ^ = 1,2 are resonantly 
coupled to the QDs. Whereas the modes with ^ 7^ 1,2 
are not within the gain spectrum of the QD ensemble or 
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have low Q-value. Since the population of the non-lasing 
modes and the cross-correlation functions 

and (&^&2) with ^^1,2 remain negligibly small, the third 
terms on the right-hand side of Eq. can be effectively 
set equal to zero. Thus, Eq. ([8]) for ^ ^ 1 , 2 can be solved 
in the adiabatic limit yielding a time constant t„; that de- 
scribes the spontaneous emission into non-lasing modes 
according to the Weisskopf-Wigner theory. The sponta- 
neous emission of QDs into non-lasing modes leading to 
a loss of excitation is described by /3-factor defined as 
the ratio of the spontaneous emission rate into the las- 
ing modes l/r; and the total spontaneous emission rate 
enhanced by the Purcell effect 1/rsp: 



/3 = 



-1 I -1 ■ 



(9) 



">P 'I ~^ 'nl 

The dynamics of carrier population of the electrons in 
the s-shell is given by 



-1 
I ■ 



(10) 

Here, the first term on the right-hand side originates from 
the interaction with the cavity-modes, the second term 
describes the relaxation of carriers from the p- to the 
s-shell with a relaxation timescale Tc and the term rep- 
resents the loss of excitation into the non-lasing modes. 

Further, we assume, that the p-shell carriers are gener- 
ated at a constant pump rate p. Then, similar to Eq. (ITOt , 
the equation of motion for the carrier population of the 
electrons in the p-shell reads: 

^(5(4cp) =p{S{vlvp) - S{clcp}) 
- S{clcp){l - d{clc,))T-'-5{clcp)il - S{vlvp))T-\ 

(11) 

where the last term on the right-hand side describes 
spontaneous recombination of p-shell carriers. The cor- 
responding equations for valence band carriers are rele- 
gated into Appendix [K\ 

The form of the expression for the intensity correla- 
tion functions suggests [see Eq. ([1])] that to exploit the 
statistical properties of the light emission using intensity 
correlations, a consistent treatment within the cluster ex- 
pansion up to the quadruplet order is required. In par- 
ticular, the equations of motion for cavity-mode intensity 
correlations read: 



dt 



+9cA^]cMblbc') + 9c'AvUM,b^)) . (12) 



The equations of motion for further correlation functions 
of the quadruplet order, which include correlation be- 
tween the photon-assisted polarization and the photon 
number, can be found in Appendix [X] [see Eqs. (|A4[) - 
^1. 



2. Results 

As described above, the quadruplet order of the clus- 
ter expansion leads to a system of coupled equations [see 
Eqs. dZ])-®, (Unil-dlll), Cai together with Eqs. SM\i - 
(|A7|) ]. The system of differential equations describe the 
dynamics of various correlations between carriers and 
cavity modes. In particular, the method makes it pos- 
sible to obtain both amplitude and intensity correlation 
functions of the cavity emission modes including the ef- 
fects of the carrier-photon correlations and the many- 
body Coulomb interaction. 

In the ensuing section the numerical analysis of the 
time evolution of the emission correlation functions is 
presented. To relate our theory to the experimental re- 
sults we estimate the number of QDs with effective gain 
contribution by starting with the initial density of present 
QDs and excluding the ones with negligible spectral and 
spatial overlap. Thus, it is assumed that the cavity mode 
field is coupled to N identical QDs. Further, we consider 
continuous carrier generation in the p-shell at a constant 
rate p as an excitation process. 

To obtain a valid comparison with the experimental 
results we simulate the coupled system using numeri- 
cal integration routines with a realistic set of param- 
eters ^ = 0.1, Ki= K2 = 0.03[l/ps], r = 2.06 [1/ps], 
Tsp — 50 [ps], Tc — 1 [ps] and t„ = 0.5 [ps]. The number 
of carriers within the frequency region of interest is esti- 
mated from the total density of QDs to be = 40. For 
the assumed /3 = 0.1 the carrier recombination is deter- 
mined by the stimulated emission into the lasing modes 1 
and 2 with a characteristic time scale t; = Tsp//3 and into 
the non-lasing modes with a characteristic time scale that 
can be found from Eq. (|9]) for the given set of parameters. 
Further, we assume that the cavity mode 1 is in exact res- 
onance with the QD transition (Ais = 0) and the mode 
2 is detuned with A12 = uji — uj2 — ^2s — 0.2 [1/ps]. In 
Fig. [3] we present the simulation results for intensity 
functions for the modes — ^ — 1,2, auto- 

correlation functions and crosscorrelation as a function 
of the pump power. Figure ^a.) reveals, that whereas 
the mode 1 shows a drastic increase of emission inten- 
sity, the intensity of the emission mode 2 reaches a max- 
imum and then slowly decreases with increasing pump 
power in good agreement with the experimental data de- 
picted in Fig. [5Ja). The calculations further show, that, 
in agreement with the experimental data in Fig. ^c), 
the dependencies of the autocorrelation functions for the 
cavity modes 1 and 2 on the pump power exhibit dra- 
matically different behavior. As shown in Fig. EJb) for 
low values of pump power, the autocorrelation function is 
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FIG. 3. (color online) Laser characteristics calculated with 
the semiconductor model, (a) Intensity correlation functions 
for the modes 1 and 2 as a function of the pump power in a 
log- log plot, (b) Autocorrelation functions of the two modes, 
(c) Crosscorrelation function between the modes 1 and 2. 



equal to 2 characteristic for the statistics of thermal light. 
For higher rates of the pump power, the autocorrelation 
function of the mode 1 drops close to the value 1 indi- 
cating the emission of coherent laser light. In contrast, 
the autocorrelation function of mode 2 slightly decreases 
at first with increasing pump powers, but for larger val- 
ues of the pump power, it increases and reaches values 
well above 2, which is in agreement with the behavior of 
the autocorrelation function detected in the experiment 
{see Fig. [Ifc); recall the limited temporal resolution of 
the HBT^''}. The gain competition behavior between 
the modes can be approved by plotting the crosscorrela- 
tion function [see Fig. ^c)], that decreases to the values 
smaller than unity at the power pump values for which 
the lasing behavior of the mode 1 is observed [also, com- 
pare to Fig. Hfd)]. Note, that the discrepancy of the ex- 
perimental and theoretical results for the autocorrelation 
function of mode 2 [Figs.[2l^c) andlS^b), correspondingly] 
and the crosscorrelation function [Figs, ^d) and EJc), 
correspondingly] at the higher pump powers is due to 
the crosstalk between the modes, which cannot be com- 
pletely avoided in the measurements. 

The numerical simulation of the cluster expansion 
truncation scheme of the quadruplet order can be ap- 
proved by plotting the emission mode autocorrelation 
functions for higher order of truncation (not shown), 



which demonstrates qualitatively the same behavior of 
the functions independent of the order of truncation. It 
is important to note, that since the framework of the mi- 
croscopic semiconductor theory presented in this section 
is based on the truncation of the hierarchy of equations 
for correlation function, the numerical results are valid in 
the regime when higher order correlations remain small. 
As it can seen from the numerical evaluation of the trun- 
cated equations, this is not the case for pump power rates 
exceeding 2 x 10~^ [1/ps], where the correlation functions 
strongly increase. To get a deeper understanding of the 
statistical properties of the emission in the next section 
we will use a different approach to gain insight into the 
full photon statistics. 



B. Extended Birth-Death Approach 

In the following subsection we present an alternative 
approach to the study of the light-matter interaction of 
QDs with a bimodal cavity, that involves numerically 
solving a complete master equation and thus deriving 
the time evolution of the system. In contrast to the 
Microscopic Semiconductor Model discussed in detail in 
Sec, nil Al the approach allows to calculate not only pho- 
ton auto- and crosscorrelation-functions but also full pho- 
ton statistics. We follow the Rice and Carmichael ap- 
proach (see Ref^i) and extend it to the case of a bimodal 
cavity with two resonant modes containing ni and n2 
photons, respectively. The method simplifies the model 
for the gain medium and takes into account only fully 
inverted two-level systems. Note, that no semiconductor 
effects or complex level structure are reflected. The state 
of the gain medium is fully described by the number of 
excited carriers N. A detailed discussion of the master 
equation approach, the semiclassical rate equations and 
its connection to the semiconductor theory for the case 
of a single-mode microcavity can be found in RefsJ^i^ 
and references therein. The master equation describes 
the time evolution of the diagonal elements 



Pn 



{ni,n2,N\p\ni,n2,N) 



(13) 



of the density matrix p. These elements can be inter- 
preted as the probability of finding a state with tii, ?i2 
photons in the modes 1 and 2, respectively, and N Atoms 
in the excited state. 



1. Physical Model 

To arrive at the final form of the master equation a 
birth and death model, analogue the one introduced by 
Rice and Carmichael^i, is considered. Transition rates 
into and out of the state p^'^^ are connected to the 
relevant processes in the coupled photon carrier system. 
Figure 0] shows how the master equation is derived on 
a phenomenological level. The filled circles represent a 
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state with TV excited carriers, rii and n2 photons in the 
cavity modes, i.e. the diagonal elements p^'"^ of the 
density matrix. The photon distribution for mode ^ is 
gained by summation over the remaining indices P{n^) — 
n( Pjv '"'"- Figure mja) illustrates the coupling of one 
mode to the gain medium. The horizontal axis shows 
the number of photons in mode ^ and the vertical 
axis shows the number of excited carriers N. The carrier 
generation is represented by solid vertical arrows since 
the photon number is not changed. The rate of car- 
rier generation in the excited level is given by the pump 
power p. Vertical dotted arrows indicate the loss of ex- 
cited carriers due to spontaneous emission into non lasing 
modes. Moreover, the emission into the cavity modes is 
represented by pairs of diagonal arrows corresponding to 
spontaneous (dotted arrow) and stimulated (solid arrow) 
emissions, since an excited carrier is lost and one pho- 
ton in one of the modes is gained. The factors Tj^^ and 
T;2^ are introduced, which represent the fractions of laser 
emission rate into the cavity modes 1 and 2, correspond- 
ingly, where the relation Ti~^^+Ti^^ = t^^ holds. Further, 
the interaction of the modes is illustrated in Fig. |H[b). 
The two axes show the number of photons ni, n2 in the 
modes 1 and 2, respectively. The horizontal and verti- 
cal dotted arrows represent the cavity losses of the two 
lasing modes with the loss rates 2k^. A phenomenolog- 
ical nonlinear coupling between the lasing modes 1 and 
2 mediated by the gain medium is also introduced. In 
contrast to the Microscopic Semiconductor Theory (see 
Sec. lIII A")) , where the couphng between the cavity modes 
is mediated by the overlap of the mode functions with the 
gain carriers, here a nonlinear coupling between the las- 
ing modes 1 and 2 is introduced phenomenologically. The 
mode coupling strengths ^12 and ^21 are represented by 
the diagonal solid arrows in the sketch and play the role 
of the detuning between the modes used in Sec. lIII Al The 
complete master equation derived by the phenomenolog- 
ical birth and death model reads: 

A "i:"2 _ r ni,n2 ni,n2] 

^^Pn —P[Pn-i Pn \ 

-r;ANpl''"-'-{N + l)pl^^n 

- r-,\n,Np'r'^ - (m - l)(iV + l)p';Vi''"1 

- r,-i[n2iVp;^'"^ - (^2 - l)(iV + l)p^^4"i^-^] 







-1; 


ni — l,n2 

Pn+1 






-1; 


ni,ri2-l 

Pn+1 




- (ni - 


1-1 


\ ni+l,n2 

iPn 




- ("2 - 


1-1 


\ ni,n2 + l 

)Pn 



- a2[nin2p^^'"^ - (ni + l)(n2 - l)p'^^+^^"^-^] 
- 6i[nin2p^^'"= - (m - l)(n2 + l)p^^-^'"^+^]. (14) 



(a) 




— 1 n( + l 



(b) 




FIG. 4. Schematic representation of the various processes in 
the extended birth-death-model. The upper sketch (a) illus- 
trates the light-matter-interaction and shows the transition 
rates into and out of the state with photon number and 
carrier number A^. Solid arrows represent stimulated emission 
and pump, dotted arrows show spontaneous emission. The 
lower sketch (b) illustrates the interaction of the modes and 
shows the transition rates into and out of the state with pho- 
ton number ni and n2. Solid arrows show mode interaction, 
dotted arrows represent cavity losses. 



In the above, each line corresponds to a process in the 
coupled carrier photon system, i.e. to arrows going in 
and out of a state p^'"'^ in Fig. S) 
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2. Results 



The stationary solution of Eq. (|T4|) gives the photon 
probabihty distribution. Figure [S] shows the photon dis- 
tributions P{n^) for various pump strengths. Above the 
laser threshold the autocorrelation functions of the modes 
1 and 2 are quite different, while ffii (0) drops to val- 

(2) 

ues close to one indicating Poissonian statistics, 522 (0) 
rises up to values substantially larger then two (thermal 
statistics). The results for the autocorrelation functions 
of the modes 1 and 2 are in full agreement with the ones 
obtained within the Microscopic Semiconductor Theory 
in Sec. IIII Al (see e.g. Fig. [3]). However, in contrast to 
the experimental data presented in Fig. [21 the autocor- 
relation function of the mode 2 monotonically increases 
also for high pump power rates. 

The full photon statistics reveals that both mode 
statistics exhibit a double peak structure. The first peak 
appears at the zero photon state declining very steep and 
a second Poissonian like peak appears at the higher pho- 
ton states. In the mode 1 the Poissonian peak is very pro- 
nounced and dominates the statistics, the mode 2 is dom- 
inated by the first peak at zero photon number states. 
These statistics combined with the fact of a crosscorrela- 
tion function far below unity allows for the interpretation 
of a switching behavior of the modes. Both modes are in 
a superposition of a lasing and a non lasing state and are 
alternating in between them. 



IV. DISCUSSION AND CONCLUSIONS 

We have investigated laser emission of electrically 
pumped quantum dots in a bimodal micropillar cavity 
with special emphasis to the effects induced by gain com- 
petition of the two orthogonally polarized modes. 

The system consisting of a single low-density layer of 
QDs and two spectrally splitted but overlapping modes 
with nearly equal Q-factors, induced from the double de- 
generate fundamental mode by slight cross-section asym- 
metry of the pillar, represents a viable platform for the 
study of the coupling of two cavity modes in the pres- 
ence of a common gain medium. The polarization re- 
solved measurements of the statistical properties of the 
emitted light reveal, that the two competing modes dis- 
play completely different features. One of the modes 
(mode 1) demonstrates typical statistical behavior of a 
laser mode, namely the mode intensity displays the usual 
"s" -shaped input-output characteristic, and the autocor- 
relation function at zero time delay, measured using a 
HBT setup, indicates the transition from spontaneous to 
stimulated emission for increasing pump rates. The mea- 
surements of the input-output characteristic of mode 2 
indicate the threshold behavior, but for further increasing 
pump rates the intensity saturates and even decreases, as 
the result of the competition of the two modes induced 
by the common gain material. Moreover, the autocorre- 
lation function at zero time delay of mode 2 at certain 



0.08 



0.06 



,0.04 



0.02 









- p=i 


gf;=2.o 








p=20 


gf,'=i.7 


1 ■ 






p =30 
















\ 
\ 




- - - p =70 


gf;=i.2 




\ 

i 




p=100 


(91 

g =1.1 
'='11 




\ 

• 


\ 










\ 

\ 








1 t 


\ , 
% / 








\ 


\ / 








% t 


/\ 








\ / 
— -V 


y \ 
1**^ * 







20 40 60 80 100 120 

"1 



(b) 



O.O81T 



0.06 



S0.04 



0.02 



I 

i 




— p=i 


g'''=2.o 


1 






^22 


i 




p=20 


g'''=2.o 








^22 


1 




p=30 


g'''=2.i 


J 






22 


n 




- - - p =70 


g,,= 3.5 


^1 








:l 

1 




p=100 


g'''=4.5 


1 : 1 






^22 


' - 1 
\l t 

































20 



40 



60 

n-2. 



80 



100 120 



FIG. 5. (color online) Photon probability distributions for 
mode 1 (a) and mode 2 (b) are shown for different pump 
rates in the units of \tJ^\ and for = = 0.05 r^"^, r~j^ 



:0.9rs/, 2kx =rsp, 2k2 = 1.2 r^/, ^12 =^21 =0.1 r^pV 



pump rates higher than the threshold values exhibits in- 
tensity fluctuation much higher than for a thermal state. 
It is worth to note, that at these rates of the injection 
current the anti-correlation between the two modes is 
the strongest. For even larger pump rates, the crosstalk 
between the modes induces a reduction of the autocorre- 
lation function at zero time delay of the mode 2 reaching 
the value for a lasing mode. Similarly, at these pump 
rates the crosscorrelation measurements indicate increas- 
ing correlation between the modes due to the crosstalk. 

The experimental results have been supported by the 
theoretical calculations within the framework of the Mi- 
croscopic Semiconductor TheoryH, which we have ex- 
tended to the case of two cavity modes interacting 
with the QD-gain medium. Using the cluster expansion 
scheme for correlation functions, we have obtained the 
emission statistics of the carrier-photon system in the bi- 
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modal cavity, taking into account the many-body effects. 
Importantly, within our approach the effects related to 
the coupling of the two modes of the bimodal cavity, in- 
duced by the interaction with the common QD carriers 
are consistently included on the microscopic level. The 
solution of the equations of motion for correlation func- 
tions reveals, that indeed the autocorrelation function of 
the mode 2 for the pump rates larger than the threshold 
rate reaches values well above ^'^^■'(0) = 2, that corre- 
sponds to the thermal state of light. The decrease of the 
crosscorrelation function of the two modes below unity 
indicates anti-correlated behavior of the mode coupling 
at these pump rates. In fact, this effect can be explained 
by random intensity switching between the two modes, 
which has negligible influence on the photon statistics of 
the lasing mode, but strongly affects the mode 2 for which 
the relative strength of fluctuations is larger. It is worth 
to mention, that in the case of macroscopic two-mode 
ring lasers22, large intensity fluctuations have been also 
found in the statistics of the more lossy mode, as the re- 
sult of the mode competition with the favored mode and 
emission switching of the common atomic ensemble. 



To complement the theoretical results with the photon- 
number statistics and to provide an interpretation to the 
super-thermal intensity fluctuations of mode 2, we have 
extended the birth-death model of the master equation 
to the case of a bimodal cavity. In particular, we have 
assumed a phenomenological nonlinear coupling between 
the cavity modes, mediated by the overlap of the modes 
with the gain carriers. In this way, solving the complete 
master equation, we have shown, that the photon num- 
ber statistics of the both modes exhibit similar double 
peak structure — a peak at the zero photon state and a 
second peak at a higher photon number. The results im- 
ply, that the both modes are in a superposition of a lasing 
state and a thermal-like state. Whereas, the Poisson peak 
at the higher photon number dominates the statistics of 
the mode 1, the statistics of the mode 2 indicates ther- 
mal state-like behavior, with a pronounced peak at zero 
photon state complemented with a local maximum at a 
higher photon number state. Thus, we may conclude 
that the photon number distribution of the two modes 
approves the switching behavior of the interaction of the 
two modes with the common gain medium. 



Similar double peak curve has been also observed in a 
semiclassical approach for the intensity probability dis- 
tributions of the both modes of ring lasers'^^, where the 
nonzero values of the light intensity are much more prob- 
able than the zero values for the favored mode and the 
other way around for the lossy mode. Moreover, a dou- 
ble peak structure of the photon number distribution has 
been found for the composite mode at threshold in the 
two-mode open laser theory^, where both modes inter- 
act with the common ensemble of atoms and with the 
common dissipation system. 
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Appendix A: Equations for Microscopic 
Semiconductor Model 



In this Appendix we present the equations of motion 
that together with Eqs. ©, dS]), (HU]), (HB) and ^ com- 
plete the full set of equations of motion for one-time cor- 
relation functions on the quadruplet level of the cluster 
expansion: 



dt 



-g^,S{clc,){l-S{vlv,))+J2 [gcAblb^)i^c]c,) ~ B,v]v,)) 
C 



^55j(5(c1.w,&^) +H.C. 



(A2) 



- (5(w^,Wp)(l - (5(wJ,w,))r,„"^-j- (5(c],Cp)(l - (5(w^Wp))r-\ 

(A3) 

- H {9i'A^^^]b\b^'bd - 9i'j^c]vjbi;)S{blb^>) 
-g^,,S(v]c,bl,blbc) - g^>,6{v]c,bl)S{bl,b^)) , (A4) 



—S{v]vjblb^) = -(kj -I- K^)5{v]vjblb^) 

-E [geMvAk'k)+9eA^hbc){i - S{v}v,) + S{blb^,)) 



+geAvhbl,blbO + geA^Ubl){l - 5{v]v,) + <5(feUc))' 



(A5) 
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-^S{c]vjblbib^^) 



c 

+ g^,jd{c]cjblb^>){l - d{v]vj) - d{bl,bc)) 

(A6) 



j^6{v]c,b\blb,) 

= i[-Ajg + i{K^ + + + iT]^v]cjb\b\b£^:) 

- gi5{c]cimv]vjb\b^,) - ^Vjblb^,) + Siblblbnb^,)) 

E [+9C'Achblb^')a - S{v}v,) + d{blbc)) 
C 

+ gcAc]c4k')i^ - ^v]vj) + S{blbc)) 
- 29cAvhbl)d{v}cjbl) - g^,jS{v}vj)d{blblbM 

(A7) 
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